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Abstract 

O 
O 

\ We compute the probabihty distribution of the invariant separation between 

nucleation centers of colhding true vacuum bubbles arising from the decay of 
a false de Sitter space vacuum. We find that even in the limit of a very small 
^ , nucleation rate {TH~^ <^ 1) the production of widely separated bubble pairs 

is suppressed. This distribution is of particular relevance for the recently 
proposed "colliding bubble braneworld" scenario, in which the value of 
(the contribution of negative spatial curvature to the cosmological density 
■ parameter) is determined by the invariant separation of the colliding bubble 

, pair. We also consider the probability of a collision with a 'third' bubble. 

o 

O : I. INTRODUCTION 

In Guth's original proposal for old inflation |[l[], a false vacuum with the geometry of de 
Ph! Sitter space decays through the nucleation as the result of quantum tunnelling of critical 
r-| ! bubbles filled with the true vacuum, presumably Minkowski space or very nearly so. The 
dynamics of false vacuum decay, which had been elucidated by a number of authors |@], may 
be characterized by the dimensionless parameter T / W^^^ (where H is the Hubble constant in 
the false vacuum and d is the total number of spatial dimensions, here d = 3). The original 
'old' inflation envisaged by Guth was unsuccessful because of the presence of two conflicting 
requirements. On the one hand, for inflation to solve the horizon problem, to smooth out 
whatever inhomogeneities may have existed initially, and to sweep away the primordial 
monopoles, T/H'^^^ must be small. Yet for inflation to end through the percolation of a 
large number of small bubbles, as Guth had originally envisag ed, r/i/'^+i had to be large. 
Although this potential conflict is noted in Guth's first paper, it was not until the more 
detailed work of Guth and Weinberg P] that this incompatibility was quantitatively and 
definitively established. 

Although this deficiency of 'old inflation' was remedied in the work of Linde 0] and 
of Albrecht and Steinhardt p[ in 'new' and 'chaotic' inflation by dispensing with bubbles 
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altogether, with inflation ending through the 'slow roll' of a scalar field along a relatively 
fiat potential, the idea of inflation with bubbles re-emerges in 'open inflation' [^, where 
re-heating occurs through the evolution of the inflaton field within the bubble interior. In 
this case T / W^^^ is so small that ony a single bubble is contemplated. 

More recently, the collision of an isolated pair of true vacuum bubbles filled with (4+1)- 
dimensional anti de Sitter space expanding within a false vacuum of (4+l)-dimensional de 
Sitter space has been proposed as a mechanism for setting up a Randall- Sundrum- like |E| 



cosmology |T0[ with well defined initial conditions A successful braneworld cosmology 
must resolve the (4+l)-dimensional smoothness problem |^ — that is, explain why the bulk 
geometry (and any additional fields in the bulk as well) respect the same three-dimensional 
spatial homogeneity and isotropy as the universe on the brane. It does not suffice to pos- 
tulate an epoch of quasi-exponential expansion on the brane itself, because after such an 
inflationary epoch ends, through their gravitational coupling, imperfections in the bulk will 
inevitably induce inhomogeneities on the brane. Because anti de Sitter space, quite in con- 
trast to de Sitter space, does not lose its hair — the amplitudes of perturbations remain 
constant rather than decay — some other mechanism to explain the spatial homogeneity and 
isotropy of the bulk, and to establish a preferred initial quantum state is required. 

In the colliding bubble braneworld scenario, an initial (4+l)-dimensional de Sitter epoch 
sweeps away whatever inhomogeneities may have previously existed. After a modest number 
of expansion times, the quantum fluctuations about the classical vacuum rapidly approach 



the Bunch-Davies vacuum [Tj]. Classically, the nucleation of a single true vacuum bubble 
breaks the 5*0(5, 1) symmetry of dS^ to >S'(9(4, 1), the subgroup of transformations that leave 
the nucleation center of the bubble invariant. For the case of two bubbles, this residual 
symmetry is further broken to S'0(3, 1), the subgroup of those transformations that 
leave both nucleation centers invariant. The three-dimensional, completely spatial surface 
of bubble collision, from which our local brane endowed with a modified (3+l)-dimensional 
FLRW cosmology arises, is endowed with the geometry of H^, of constant negative spatial 
curvature. 

Owing to the negative spatial curvature of this geometry, the eventual fate of such a 
universe (in the absence of a non-vanishing cosmological constant or something similar) is 
that of an empty expanding universe whose expansion is predominantly driven by its spatial 
curvature, as in a Milne universe. However, because of the exponential suppression of the 
nucleation of bubbles, this dreary fate may lie far in our future. It is not difficult to envisage 
bubble separations so large that today the spatial curvature is still negligible. f2o, the 
cosmological density parameter today (when T^ad = '2.75K), is a monotonically increasing 
function of the bubble pair separation. The farther apart the bubbles, the larger pcou and 
the smaller the spatial curvature on the surface of collision. The radius of spatial curvature 
of the surface of collision is given by |]13| 



R = Hds-'tan[a/2], (1) 

which diverges as a ^ vr. Here cr is the spatial geodesic separation between the nucleation 
centers in units where H^^ = 1, so that a = n corresponds to bubble separated exactly 
at the threshold of never striking each other. We observe that even though cr for colliding 
bubbles is bounded from above by a finite limit, the surface of collision can be arbitrarily 
close to flat. The gamma factor 7 = 1/^1 — (f/c)^ of the colliding bubbles in the local 
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center-of-mass frame of the collision, which is equal to 



_ sin[(T/2] 
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where is the radius of the critical bubble (in units in which H^g^^ = 1), however, is 
limited by a finite upper bound.Q Both effects drive Qq toward unity as a ^ vr. The precise 
relationship between a and Qq depends on the details of the equation of state after the 
collision, which is unknown and model dependent. 

In the colliding bubble scenario the value of Q today is a random variable dependent on 
the interbubble separation which is stochastically determined. The object of this article is 
to study the distribution of the interbubble separations for the case of bubbles nucleating 
in de Sitter space. 



II. BUBBLE COLLISION PROBABILITIES 

We now compute the probability distribution for the interbubble separation of pairs of 
colliding bubbles. This separation a is defined in a coordinate free manner as the length 
of the (shortest) spacelike geodesic connecting the nucleation centers of the two colliding 
bubbles. We adopt units in which H = 1 and assume a spatially and temporally uniform 
nucleation rate F, the absence of correlations of the nucleations of nearby bubbles, and that 
bubbles cannot nucleate inside regions already within bubbles. 

First a lightning review of the geometry of de Sitter space. Maximally extended (d+1)- 
dimensional de Sitter space dS^'^'^^^ is most readily constructed as the embedding of the unit 
hyperboloid defined by 

- + Xo' + . . . + X/ = 1 (3) 

in (d+2)-dimensional Minkowski space M^'^'^'^^ with the line element 

ds'^ = -dT^ + dXo'^ + ... + dX/. (4) 

This space in its entirety may be covered by the coordinates 

T = sinh[r], 
Xq = cosh[r] uq, 

Xd = cosh[r] rid, (5) 

where (no, . . . , rid) are the coordinates of a point on the unit sphere S'^ embedded in E^'^'^^\ 
with the line element 

ds'^ = -dT'^ + cosh^ [r] dnfa) . (6) 



^Note that the notation here differs from that of ref. IT 
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Any two points with a spacehke separation and the property that their forward hghtcones 
eventually intersect may be transformed to lie on the throat of the hyperboloid defined in 
eqn. (i.e., the sphere T = 0). For the largest such spatial separation, a = n, and two 
points are antipodal. This situation (with our idealization of zero critical bubble radius) 
corresponds to bubbles which almost but never collide. Larger spatial separations are also 
possible, but in this case no spacelike geodesies exist linking the two points. 

As a technical device, it is useful to characterize the invariant separation between two 
points P and Q in terms of the (d+2)-dimensional Minkowski space inner product 



/(P, Q) = -TpTq + J2iX^)p (X,)q. 



i=0 



It follows that 

/(P, Q) = cos[cr] 
for spacelike separated points of the first kind and that 

/(P,Q)<-1 

for those of the latter kind. (Timelike separations correspond to /(P, Q) > +1.) 



(7) 



(8) 



(9) 




FIG. 1. The spatially flat coordinate slicing covers precisely half of maximally extended 
de Sitter space, as indicated in the Penrose- Carter conformal diagram above. The two vertical 
boundaries are identified and on the horizontal boundaries the conformal factor diverges. 

In terms of the fiat coordinates, which cover precisely half of maximally extended de 
Sitter space, as indicated in Fig. 1, the line element is 



ds^ = -dr + e 



2t 



(10) 



dxi + . . . + dxd 

which after the change of variable to conformal time r] = — e^*, with — oo < r] < 0, becomes 

1 



ds' 



7]^ 



—dr] + dxi + . . . + dxa 



(11) 



The embedding of these coordinates into M^'^^^^ is given by 
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T = sinh[t] + e* 
Xq = cosh[t] - e* xV2, 

X,; = 



6 Xi 



where x ■ y 



xiyi + . . . + Xdiid- The invariant separation defined above is given by 



(12) 



(13) 



We now turn to the issue of initial conditions. It is necessary to specify an initial surface 
over which there is everywhere false vacuum. As has been previously pointed out in the 
literature 
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an unstable de Sitter vacuum extended infinitely far into the past of maximally 
extended de Sitter space does not make any sense. Bubbles formed during the contracting 
phase eventually percolate. However, if consideration is restricted to a subregion of de Sitter 
space over which there is only expansion (such as that covered by the fiat coordinates), then, 
if r is sufficiently small, infiation is eternal in the forward time direction (i.e., the bubbles 
never percolate), and presumably any conclusions drawn are insensitive to the character of 
this initial surface. 

We employ the fiat coordinates, assuming a false vacuum everywhere on the surface t = t_ 
and considering only bubbles that nucleate for t in the interval t_ < t < t^. (Eventually we 
take the double limit t_ — oo, +oo.) 




FIG. 2. We take as given that a bubble has nucleated at the point P. The surface S indicates 
the surface of possible nucleation centers separated by the same proper distance a from P. For a 
bubble that nucleates at Q on S to collide with the bubble that nucleated at P (rather than with 
some other intermediate bubble), it is necessary that no bubble nucleated in the singly shaded 
region. 

The probability density for a bubble to nucleate at time t among all the bubble nucleating 
in the time interval t_ < t < t+ is given by the normalized density 



dPf 



bn 



dt 



it) 



dt' e"'' e-^^-(*') 



(14) 



where Vic is the volume of the past light cone of the nucleation center lying above the initial 
surface t = Here we have considered an infinitesimally thin tube of constant co-moving 
d- volume, whose physical d- volume scales as e*^*. 
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Next we consider bubble collisions. At first we restrict ourselves to the simplest case of 
bubble in (l+l)-dimensional de Sitter space to avoid the complications arising from the fact 
that in higher dimensions a bubble can collide with many other bubbles. However, for all 
volumes and measures we retain the general formulae, valid for all d. In (l+l)-dimensions, 
each bubble collides with precisely one bubble from the left and one bubble from the right. 
Once the (l+l)-dimensional case has been thoroughly examined, we generalize to arbitrary 
dimension in the small F limit. 

We now calculate the conditional probability dPbc{<y\t) / da that, given that a first bubble 
has nucleated at time t, another bubble whose nucleation center is separated from that of the 
first bubble by a proper distance a collides with the first bubble. This probability density 
is given by the expression 

= / rfS r e-r^-'(^). (15) 



da Jt. 

The various volumes are indicated in Fig. 2. P is the nucleation center of the bubble that 
has by assumption nucleated at time t. The presence of this bubble implies that no other 
bubble has nucleated in the part of its past lightcone to the future of the initial surface (i.e., 
the doubly shaded region of the figure). S is the locus of all points a proper distance a to 
the right of P and rfS is the volume element on this surface. Q is a point on this surface, a 
potential nucleation center from a bubble a proper distance a from P that may collide with 
the bubble emanating from P. However, for such a collision to occur, it is necessary that no 
bubble nucleate in the lightly shaded region, of volume Vexd- It is necessary to exclude the 
case of a bubble nucleating in the volume Vexd, for any such bubble would collide first with 
the bubble emanating from P. The probability of no nucleation in the lightly shaded region 
is exp{—TVexd), hence the appearance of this factor in eqn. ([T5|) . If we remove the upper 
cut-off, taking t+ — > +00, the integral of the density from a = to a = vr becomes equal to 
unity — in other words, a collision occurs with unit probability. 

We may combine eqns. (|T^ and ([1^) to obtain the probability density properly averaged 



over the nucleation time for the first bubble t for collisions of bubbles separated by a proper 
distance a 

dPav{(r) ,dPbc, ... dPbnit) 

dt -7-{a\t) —J—. (16) 



da Jt^ da dt 

Particularizing to (1+1) dimensions and replacing the proper time t with the conformal 
time rj, we rewrite eqn. ([I5| ) as 



dP r r ^ /"'rfx(5(cos[(T] -cos[(t]) e-^^-'. (17) 



(icOs[(T] J-00 1]"^ Jo 

Without loss of generality we may set the coordinates of P to r/p = —l,xp = 0. We eliminate 
the integration over x by considering 

. rn T]^ + r]p'^ - {x - X pf - 7/^ - 1 

fix) = COS a = = , (18) 

2r] rjp ZT] 

so that eqn. (|T7|) becomes 
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dr] 



df 



dx 



^ ro dri /-r/\ ^_rK,4.(,), ^) 
J-oo -rf \ X J 



V 

dr] 



-rVe^ci(a;(r;), r;) 



^1 Jl + 2r]cos[a]+r]^ 



(19) 



We now turn to the computation of Vexd{x,ri). As can readily be seen from Fig. 2, this 
volume depends only on the point of intersection (x, f/) of the left-moving null geodesic 
through (x, 7]) with the right-moving null geodesic emanating from P. It follows that 



1 r 



^=2 



1] + X 



We compute 

Vexcliv) 



"§2(1 -1^1) + 



-'"^2(-|^|-„) = 2hJi 



1 r]^ 



(20) 



(21) 



Consequently, eqn. ([T9| ) becomes 



rfcosM 



r 



d7] 



1 



^1 Jl + 2r7cos[a] +r/2 



exp 



+2r log 



1 — f] — Jl + 2r] cos[(j] + 1] 



+ 00 



1 



^7 ^Jl — 21] cos[cr] + r]^ 



l+r/— -y/l — 2?] cos[cr] + r]^ 



2r 



(22) 



For small F, the dominant contribution to the integral lies in the neighborhood of ?7 = 0. 
Therefore, we may approximate the above expression as 



F 



F 



+ 00 



dr] 



1 



V wl — 2?] cos[cr] + r]^ . 



r]{l + cos[o-]) 



1 2r 



1 + cos [cr] 



1 /I + cos[o-] 

2 \ 2 



2r 



2r 



dr] 7] 
1 

2' 



(2T-1) 



(23) 



In other words, for F 0+, we obtain a uniform distribution in coso", or 

dP sin [a] 
la ~ 2 



(24) 



We now generalize to d > 1. Since any given bubble suffers an indeterminate number 
of collisions with other bubbles, it is not at the outset clear how to weight the various 
collisions to obtain a sensible distribution for dP/da. We first define the distribution in a 
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rather artificial way, later showing that the result obtained is more generally valid in the 
small r limit. 

We start by postulating as given the nucleation of a bubble, taken without loss of gen- 
erahty to be situated at 77 = — 1, x = 0, and compute the distribution for the first other 
bubble to coUide with this bubble. (This may seem a bit artificial because a preferred time 
direction has been singled out.) Defined in this way, eqn. (19) becomes 



Vexclilvl) 



d-i 



d 



+ 



1 dr) 

M Vd+i 

°° di] 



Vd+i 



{r) + l-2\fi\)''-{l-v) 
\rj + l-2\r]\f-{r]-ir 



(25) 



where fld-i is the area of the (d-l)-unit sphere. In the small T limit, we are only interested 
in the most singular part of Ve^ci as I77I — > + . In the limit, the second integral in the 
expression above is finite, and consequently can be ignored. The first integral is dominated 
by the contribution in the neighborhood of the lower limit of integration, which may be 
approximated as 



^excl ~ ^ 

We obtain (setting F = 



drj 



dP{a) 
dcos[a] 



(d-i)) 
Jo 



dr] 



d+l 



2(77 -|?7|) 



-d+l 



d{d-l) 



(26) 



: exp 



dr] 

r / ^ exp 
Jo ?7" 



^1 — 2r]cos[a] + r] 

2f f r]{l + cos[a]) 



i' r+00 

r / dz exp 

— 1) Jo 





d{d-l] 

d+l" 
-d+l 



d{d-l) 



d-1 





'a' 




^COS^ 




) 




.2. 





so that 



d / 1 + cos[cr] 

— — - sm cr 

da 2 \ 2 



d-1 



(27) 



dP{a) 



d-1 



(28) 



We observe a power low suppression for large bubbles. 

We now explain why the distribution above is more broadly valid, when multiple collisions 
have been taken into account. As before, we consider a bubble whose nucleation center has 
been placed at = — 1, x = 0. Wc consider an infinitesimal region on the surface of the 
bubble from (f^, 1 + f/) and (f/ + df], 1 + f] + dff) and some interval of solid angle as well, and 
consider, given that some exterior bubble first strikes this bubble in that patch, what is the 
distribution of a. Computing this distribution involves extending the patch radially outward 
and into the past along the relevant null geodesies. Let ^ range from to +00. Then 
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r = 1 + f] + ^. 



(29) 



It follows that 



cos cr 



27] 



(30) 



Therefore, ^ = —fjta.n'^[a/2]. For small F, we may take f] to be very close to zero, because 
all but a few bubbles nucleate very long after the nucleation of the first bubble. Using the 
fact that 



dP 



dV 



dcos[a] dcos[a] 



we obtain 



dV ^d^ 



d+1 



di_ 

.d+1 



di 



7]' 



d+1 





'a' 




^COS^ 




) 




.2. 





d{cos[a]) 



(31) 



(32) 



confirming the result in eqn. (^). By means of this last calculation, all bubble collisions 
for which a portion of the bubble collision surface is tangent to one of the t = {constant) 
surfaces are counted and given equal weight. 



III. COLLISIONS WITH A THIRD BUBBLE 

In the colliding bubble braneworld scenario the universe that arises from the collision of 
two bubbles has a hyperbolic geometry that would be infinite in extent in the absence of 
any collision with a third bubble. However, it is inevitable that sooner or later collisions 
with third bubbles occur, and such collisions truncate the extent of this hyperbolic universe. 
Physically, the relevant question is whether such a cosmological scenario is likely to contain 
a hyperbolic patch large enough to contain the entire universe observable to us today, for 
the consequences of a collision within our past lightcone with a third bubble would hardly 
be subtle. At a very minimum, such a collision would create an 0(1) perturbation in AT/T 
of the CMB anisotropy. More dire consequences are also possible. 

Concretely, to obtain an estimate of the likelihood of a collision with a third bubble 
within our past lightcone, we consider how the infinite surface of collision of the 

two bubbles is truncated by collisions with other bubbles. We take as given the initial 
collision of the two bubbles, assuming at least that the point ^ = on H^'^"^^ has been 
spared from a colhsion with a third bubble. Here ds'^ = B^{d^'^ + sinh^[^](if2^^_^^^ is the 

metric on H^'^~^\) We calculate the probability distribution for the radius ^ of the largest 
sphere = {(^, 6',0)|^ < ^} C H^'^~^^ such that no bubble has nucleated within the past 
lightcone of B^. As before we employ the coordinates ds"^ = [1 / 7f){—d7f + dy?) for the larger 
{d + l)-dimensional de Sitter space into which H^'^~^'> is embedded. If the integration over 
the volume contained within a past lightcone of a spacetime point is taken all the way back 
to T] = — oo, the resulting integral diverges logarithmically as ln^|?7mj„|^. However, since by 
assumption no bubble has nucleated in the past lightcone of the origin O at ^ = 0, we are 
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interested in the volume that hes between the past hghtcone of and the past hghtcone 
of the origin O, and the integral defining is convergent. One obtains 



dP dV r 

di di L 



(33) 



For a typical value of ^, TV{^) ~ 1. 

We first consider the simplest case with (c? = 2), later straightforwardly generalizing to 
d > 2. We define the hyperboloid H^'^~^^ as the locus of points a proper time r to the future 
of its focus, where proper time is calculated using the de Sitter metric and sinh[r] = R. 
Without loss of generality we fix the focus to the point t] = —1, x = 0. From the relation for 
the invariant separation in de Sitter space, it follows that the hyperboloid is defined by the 
expression 



-2r] 



-2f] 



cosh[r] = Vl + -R^. 



(34) 



Here [r], 0) is the spacetime position vector of the origin of the hyperboloid O. It follows 
that -R = (1 — ff )/{—2ff). We find it convenient to parameterize the hyperboloid according 
to 



m 

r{e) 



Rtaxi[e], 



where < 6* < Omaxi such that tan[6'm,ax] = 1/-R. As i? — > oo, O^ax 0. 

In terms of the distance in the larger de Sitter space, the line element along the 
boloid gives 



ds 



Rsec[9] 



de Vl + R^ - Rsec[9y 
or in terms of the natural dimensionless distance on H^'^~^^ 



di 

dd 



^1 + R^ - Rsec[e] 

To compute dV/dC, in eqn. (|33|) at a certain value of C,, we use the relation 

dV _ de dV 

~di ~ d^He 

and compute the infinitesimal volume {dV/d9)d9 swept out by the displacement d9 
past hghtcone at 9 is displaced to {9 + d9). In the underlying Minkowski space, with 
ds'^ = —drf + rfx^, the cone at P{9) = {ri{9),r{9),0) suffers a rigid displacement by 



_d_ 
89 



Rsec[9] ■ 



tan[^]- + sec[^]- 



We parameterize the surface of the cone using w, Wmin < w < +oo, and ip, —tc < ip 
where Wmin = -'7(6'), 



(35) 
hyper- 

(36) 
(37) 

(38) 

as the 
metric 

(39) 

< +vr. 
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V 

r 

X 



[w 



(w — w„ 



COS 



sm 



(40) 



The volume swept out by a surface element {dw){d%p) is given by the determinant 





dV 



dtpdw dO 



-Rsec[6'](w — 



taxi[e] sec[^] 
—1 cos [■?/'] sm.[ip\ 
— sin[?/'] cos ['?/'] 

R?,ec[6]{w — Wrnin) tan[^] + sec[6'] cos| 
_Rsec^[^^](w — Wmin) COS ['?/'] — cos[7r/2 + 9] 



(41) 



Since we are interested only in the volume swept out in the forward direction, we restrict 
the integration to the subregion over which this determinant is positive — in other words, 
over the range < 7r/2 + 6. Including the conformal factor to reflect de Sitter rather 
than Minkowski volume, we obtain 



dV 

'de 



2Rsec^[9] 



dw 



Me)\ w 



VT+W - Rsec[9] 



cos [6*] + 



n/2+e 



dip COS ['?/'] — cos[7r/2 + 9] 



TC 



+ 9] sm[9] 



(42) 



which in light of eqn. (pT]) becomes 

-R 



- + 9)tan[9] 



(43) 



We are concerned with the limit F <^ 1, in which the values of 9 of interest lie very near 
9max- In this limit, r 1 and tan^^ ^ 1/R. For large R, 9max is small, and dV/dC, ~ R, 
giving V{^) ~ R^. Since F is exponentially small, the requirement that TV{{^)) = 0(1), 
where here (^) is the expectation value of ^, implies that (^) is exponentially large. 

To generalize to higher dimensions, we must modify eqn. (^) by changing the conformal 
factor from to w~^^~^^'^ and replacing the factor of = 2, with where f2i = 2n, 

Q2 = 47r, .... Moreover, dr becomes r'^~^dr, where now the position of the origin for r does 
matter. With these modifications, eqn. (^21) is now 



^ = fi(„_2)i?sec2[0] P -%iw - m\) 
d9 J\vW\ 



TT/2+e 



di/j cos[V^] - cos[7r/2 + 9] r{9) + (w - rj{9) 



cos 



d-2 



^](rf_2)^^sec2[e] 



+°° du 



d+l 



(u-l) 



u 



TT/2+e 



dip cos ['?/'] — cos[7r/2 + 9] 



r{9) 
W)\ 



-I d-2 



+ (u — 1) cos[-?/'] 



(44) 



where u = w/\ri{ri)\. 
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For the dimension of greatest interest (d = 4), it is not possible to evaluate 

dV 



Vie) 



ad 



(45) 



with the integrand defined in eqn. (44). However, for F ^ 1, we are interested in 6 very 
near 6max, so that r{9) /\ri{6) \ ^ 1. Hence, we may drop all but the leading term in powers 
of r{6) /\ri{6)\ in eqn. (0), which will then reduce to 



cos[9] 



6 sm[e] 



r — > 1 and 



so that 



V 



TT 



R+ 2 +^-ax 



From eqn. (|37| ) it follows that 



1 cos[^] sin[6'm 

\r]{e)\ ~ cos[e] - cos [ft 

cos [^rnaa 



cos[6'^ 



d9- 



'max 
1 



o 



Therefore, 



and 



d9 



^ ^ -ln[6'„ax - 0]. 

(Omax -9) ^ exp [-i] 



and 



V 



TT 



-(-R + — + Omax) 



i?2 



exp [+2^] 



Given that [|], [|T3| 



F ^ exp[-y4(m4£)2] 
where A = 0(1), we find that for {m^i) ^ 1, 



(46) 



(47) 



(48) 

(49) 

(50) 
(51) 

(52) 
(53) 
(54) 



where we have suppressed all logarithmic corrections. This result implies that for large 
{rriit), the 50(3,1) symmetric regions where the bubbles collide are composed of patches 
of a size almost always containing several curvature lengths, suggests that one is unlikely to 
observe a collision with a third bubble. 
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IV. DISCUSSION 



The object of this study was twofold: first, to determine whether in the hmit F — >■ 
it is possible to obtain, with a non-negligible probability, bubble separations very near but 
just short of the maximal separation a = tt, where the bubbles just barely strike each other; 
and, second, to determine whether collisions with a third bubble constrain or rule out the 
colliding bubble braneworld scenario. 

With respect to the first question, we found that as F — > 0+, the probability distribution 
dP/da approaches a limit independent of F, given in eqn. (^), having a power law suppres- 
sion of values of a near tt. Consequently, no matter how small F, R (the curvature radius of 
the surface of collision H^) is typically of order H^g, the curvature radius of the de Sitter 
space into which the bubble expands. To obtain a value of Qq close enough to one today, 
H^g must be sufficiently large. As discussed in ref. fl^, estimating Rmin depends on the 



unknown details of the equation of state on the brane arising from the collision, although in 
a rather insensitive way. For purposes of illustration, we adopt the assumption of a radiation 
equation of state in the aftermath of the collision, obtaining^ R ^ 10^^ i {m4i)~^'^^^\ Here 
£ is the AdS curvature radius inside the bubble. The largest admissible value i ^ 1 mm 
gives R > 10^° cm. While this large curvature is somewhat of an embarrassment, the five- 
dimensional bulk cosmological constant problem implied here is orders of magnitude milder 
than the four-dimensional effective cosmological constant today. One might hope to evade 
this restriction by appealing to the anthropic principle. In other words, one would have H^g 
small so that most bubble collisions would result in empty and uninhabitable universes. We, 
however, would descend from the collision of one of those rare pairs for which Qo is not 
minuscule. Such a scenario, however, would predict a distribution of Qq seen by us peaked 
near the lowest acceptable value, much as in Weinberg's proposed anthropic resolution of 
the cosmological constant problem [^, and this prediction is at variance with observation. 

With respect to the second question, we found that collisions with third bubbles are rare, 
allowing the surface of collision to span several curvature radii, a size in excess of that of 
the universe observable to us today. Our analysis, however, considered only bubbles that 
collide with the surface of collision. Bubbles that collide from the side with one of the two 
bubbles from which our universe originated do lie in our past lightcone and in principe could 
affect us. In the idealised colliding bubble scenario such a collision simply results in another 
universe like ours; however, the possibility that some of the debris from the collision escapes 
into the bubble interior and subsequently propagates to our universe cannot be ruled out. 
We have not considered such possible collisions. 
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Note that the directions of the inequalities in eqns. (66), (69), and (70) of ref. |13| should be 
reversed. 
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